Acessibilidade / Reportar erro

Prediction of Equipment Effectiveness using Hybrid Moving Average-Adaptive Neuro Fuzzy Inference System (MA-ANFIS) for decision support in Bus Body Building Industry

Abstract

Managers are driven to accomplish significantly higher levels of operational performance due to the difficulty of today’s dynamic production environment. Typically, the precision of production facilities and the efficiency of manufacturing systems are significant variables in productivity. Thus, predicting machine performance has become an inevitable challenge for production managers. However, the question of how managers can reliably assess the effectiveness of equipments for resource allocation remains unaddressed properly. This issue has received little attention in previous research, but it is important in today’s manufacturing environment. This study introduces a hybrid moving average - adaptive neuro-fuzzy inference system (MA-ANFIS) to predict the possible effectiveness of equipment. Three real-world problems are considered when developing and evaluating three distinct equipment effectiveness prediction models. The evaluation confirms that the hybrid MA-ANFIS model based on Gaussian membership function outperforms other developed models. This comprehensive solution is packaged as a decision support system. This aids production managers in evaluating the equipment effectiveness, and effectively improving equipment’s performance to reduce time and cost of bus body building.

Key words
adaptive neuro-fuzzy inference system (ANFIS); moving average (MA) model; hybrid extremal-micro genetic algorithm (Ex-µGA); predictive analysis; equipment effectiveness

Introduction

The competitive industrial environment forces industries to maintain consistent productivity. Therefore, maintenance activities are becoming increasingly important to enhance equipment performance and ensure consistent production. Total Productive Maintenance (TPM), Just In Time Manufacturing (JITM) and Business Process Reengineering (BPR) are just a few technological examples of machine performance improvement in the industrial sector (Ljungberg 1998LJUNGBERG Õ. 1988. Measurement of overall equipment effectiveness as a basis for TPM activities. Int J Oper & Prod Manag 18(37): 510-528. URL http://dx.doi.org/10.1108/01443579810206334%0., Nikolopoulos et al. 2003NIKOLOPOULOS K, METAXIOTIS K, LEKATIS N & ASSIMAKOPOULOS V. 2003. Integrating industrial maintenance strategy into ERP. Indust Manag Data Syst 103(3-4): 184-191. doi:10.1108/02635570310465661., Chou et al. 2005CHOU TC, HSU LL, YEH YJ & HO CT. 2005. Towards a framework of the performance evaluation of SMEs’ industry portals. Ind Manag Data Syst 105(4): 527-544. doi:10.1108/02635570510592406.). The TPM is widely considered as the most promising strategic technology for assisting industries in achieving competitiveness (Nakajima 1988NAKAJIMA S. 1988. Introduction to total productive maintenance. Student Study Guide. New York: Productivity Press, p. 129.). In practice, TPM is an integrated approach based on American and Japanese maintenance principles for maintaining equipment performance and reducing machine failures. This also includes employee participation in self-maintenance (Ahuja & Khamba 2008AHUJA IPS & KHAMBA JS. 2008. Total productive maintenance: literature review and directions. Int J Qual Reliab Manag 25(7): 709-756.). The quantitative metric assessment of the TPM process named as Overall Equipment Effectiveness (OEE) (Anvari et al. 2011ANVARI F, EDWARDS R & STARR A. 2011. Evaluation of overall equipment effectiveness based on market. J Qual Maint Eng 16(3): 256-270.). The OEE is an issue that draws the attention of production managers for the purpose of enhancing productivity. In this regard, continuous OEE monitoring is carried out in order to assess equipment efficiency and propose appropriate measures for its management. Benchmarking analysis is commonly used to determine the best management techniques for increasing manufacturing plant efficiency and the effective rate of OEE (Fattah et al. 2017FATTAH J, EZZINE L & LACHHAB A. 2017. Evaluating the performance of a production line by the overall equipment effectiveness: An approach based on best maintenance practices. Int J Eng Res Afr 30: 181-189. doi:10.4028/www.scientific.net/JERA.30.181.). In particular, the OEE is quantified by combining three distinct variables into a single one. The variables are equipment availability (A), equipment performance (P), and quality (Q) of production (Sharma 2019SHARMA R. 2019. Overall equipment effectiveness measurement of TPM manager model machines in flexible manufacturing environment: a case study of automobile sector. Int J Prod Qual Manag 26(2): 206-222. doi:10.1504/ijpqm.2019.10018982.).

There have been significant progresses in assessing OEE, since then Wang & Lee (2001)WANG FK & LEE W. 2001. Learning curve analysis in total productive maintenance. Omega 29(6): 491-499. doi:10.1016/S0305-0483(01)00039-1. formulated time constant regression model. In that study, expected OEE is used to compare the maintenance performance of TPM implementations that have been implemented in two different industries from Taiwan and Japan. The findings imply that both industries perform similarly in terms of maintenance, and that expected OEE may be easily determined by tracking maintenance progress. However, the equipment uncertainties are not taken into account in the study’s input parameters. Ma et al. (2012)MA YH, ZHENG YB, ZHANG Y, JIANG Y & ZHAO LL. 2012. Diagnosis and prediction for equipment efficiency based on OEE in FCC M6 welding workshop. Appl Mech Mat 130-134: 1092-1096. doi:10.4028/www.scientific.net/AMM.130-134.1092. used exponential smoothing technique to measure the OEE for a welding workshop. The results confirmed that the smoothing time series observation has the best effect while considering loss functions of A, P, and Q. Relkar & Nandurkar (2012)RELKAR AS & NANDURKAR KN. 2012. Optimizing & analysing overall equipment effectiveness (OEE) through design of experiments (DOE). Procedia Eng 38: 2973-2980. doi:10.1016/j.proeng.2012.06.347. developed a regression equation for predicting OEE using Design of experiments (DoE). The study found that focusing on equipment performance would improve OEE. Likewise, Yuniawan et al. (2014)YUNIAWAN D, ITO T & MOHAMAD ED. 2013. Measurement of Production Line Performance Based on OEE Extension. In: The Proceedings of Manufacturing Systems Division Conference, 57-58 p. doi:10.1299/jsmemsd.2014.57. integrated the Arena simulation model with DoE technique. The capability of the suggested anticipated model enables decision makers to benefit. Regrettably, the two inquiries do not provide a systematic approach for defining the input parameters. In addition, Kurscheidt Netto et al. (2017)KURSCHEIDT NETTO RJ, SANTOS EAP, DE FREITAS ROCHA LOURES E & PIEREZAN R. 2017. Using Overall Equipment Effectiveness (OEE) to Predict Shutdown Maintenance. In: AMORIM M, FERREIRA C, VIEIRA JR M & PRADO C (Eds), Engineering Systems and Networks. Lecture Notes in Management and Industrial Engineering. Berlin: Springer, p. 13-21. doi:10.1007/978-3-319-45748-2_2. used a statistical process control tool to link cycle time and equipment stoppage trends. This finding demonstrates that equipment breakdowns are putting a strain on the equipment’s performance, resulting in an inaccurate prediction of OEE. Nonetheless, this study has focused solely on the performance loss function. In a another research Prasetyo & Veroya (2020)PRASETYO YT & VEROYA FC. 2020. An Application of Overall Equipment Effectiveness (OEE) for Minimizing the Bottleneck Process in Semiconductor Industry. In: 2020 IEEE 7th International Conference on Industrial Engineering and Applications, ICIEA 2020, p. 345-349. doi:10.1109/ICIEA49774.2020.9101925. investigated continuous improvement tools to identify bottlenecks in the semiconductor sector in the Philippines in order to improve OEE. The results show a 30 percentage increase in OEE and recommend employing the DMAIC method to improve OEE in the long term. The main problem with these models is that there is no systematic method for determining the membership function parameters, which must be determined by expert knowledge of the represented system.

Furthermore, substantial progress has been made in analyzing OEE using machine learning. In the context of machine learning, Ok & Sinha (2006)OK SC & SINHA SK. 2006. Construction equipment productivity estimation using artificial neural network model. Construction Management and Economics 24(10): 1029-1044. doi:10.1080/01446190600851033. developed a Neural Network model for predicting construction equipment performance. In this study, seven different factors were taken into account. The results were compared to the results of a multiple regression model. The Neural Network (NN) model outperforms the others with a minimal mean square error of 0.0014. Likewise, Kuo & Lin (2010)KUO Y & LIN KP. 2010. Using neural network and decision tree for machine reliability prediction. Int J Adv Manuf Technol 50(9-12): 1243-1251. doi:10.1007/s00170-010-2593-7. developed four different Neural Network models for equipment performance forecasting in a washing machine manufacturing company. According to the comparison analysis, the combination of NN and decision tree model is the most likely to reproduce reality, with an average test data error of 8.73 %. The NN is a set of algorithms meant to recognize patterns that are similar to those seen in the human brain. These can ability to learn from a set of input and output data in real time. The adaptive neuro-fuzzy inference system (ANFIS) is gaining significance in dealing with inadequately and unexpected domains such as OEE forecasts as a result of these insights. In a research Bekar et al. (2015)BEKAR ET, CAKMAKCI M & KAHRAMAN C. 2015. An ANFIS algorithm for forecasting overall equipment effectiveness parameter in total productive maintenance. J Mult-Valued Log Soft Comput 25(6): 535-554. used the concept of ANFIS for prediction of OEE in production system. However, triangular membership function was solely used for determining of the membership degrees for each input parameter. Other membership function was not taking in account for finding membership degrees. In addition, Hassani et al. (2019)HASSANI IE, MAZGUALDI CE & MASROUR T. 2019. Artificial Intelligence and Machine Learning to Predict and Improve Efficiency in Manufacturing Industry. arXiv URL http://arxiv.org/abs/1901.02256.
http://arxiv.org/abs/1901.02256 ...
proposed five different machine learning (ML) algorithms for forecasting expected OEE values. According to the findings, deep neural networks perform much better, with a mean absolute error of 6.27 and a mean absolute percentage error of 11.76%. However, this prediction model is only a test case. A combination of genetic algorithm optimization must be incorporated to build a prediction model for real-time case validation. Likewise, Engelmann et al. (2020)ENGELMANN B, SCHMITT S, MILLER E, BRÄUTIGAM V & SCHMITT J. 2020. Advances in machine learning detecting changeover processes in cyber physical production systems. J Manuf Mat Process 4(4): 108. doi:10.3390/jmmp4040108. proposed a machine learning model and evaluated it using the MATLAB classification learner to forecast OEE based on machine availability. Based on the results, the fine tree decision algorithm provided the greatest fit, with an overall accuracy of 92.8%. The RUS boosted tree, on the other hand, outperforms the fine tree by 0.93 to 0.92. As a result, the RUS boosted tree technique was found to be superior in classifying the changeover process. The model provided here solely considers machine availability losses. These findings indicate that when developing OEE prediction models, a more thorough assessment of input parameter selection is required. However, the approach is determined by the type and volume of available input data.

According to the literature, ANFIS has been recognized an improved prediction tool than other methods. The learning abilities of NN were combined with the reasoning abilities of fuzzy logic to build a working platform. ANFIS provides specialized knowledge in the form of fuzzy "if-then" rules based on a membership function approximation from specified determinants and response data sets. The fuzzy set theory handles the uncertainties associated with network computation, while the NN provides model adaptability. The membership function parameters linguistic rules can adjust directly from NN training capabilities concerning refining model performance. The development of a hybrid supervised learning system that can predict the relationship between predictor and responder variables using artificial neural abilities based on evidence-based intelligence. Therefore, the hybrid ANFIS approach is selected to propose an OEE prediction model.

The computation of OEE allows managers to keep track of their processes and pinpoint the major losses that limit machine effectiveness using substitution of latest technology like IoT. However, this estimate is posses huge cost. The goal of this work is to develop a hybrid ANFIS model that can forecast the estimated OEE value, allowing managers to assess the equipment’s effectiveness with low cost. In this study, 3 input parameters (production losses) are selected to develop prediction model from the bus body building industry. The moving average smoothing technique is used to smooth the input data based on the findings of Ma et al. (2012)MA YH, ZHENG YB, ZHANG Y, JIANG Y & ZHAO LL. 2012. Diagnosis and prediction for equipment efficiency based on OEE in FCC M6 welding workshop. Appl Mech Mat 130-134: 1092-1096. doi:10.4028/www.scientific.net/AMM.130-134.1092. and Lotfi et al. (2020)LOTFI K, BONAKDARI H, EBTEHAJ I, DELATOLLA R, ZINATIZADEH AA & GHARABAGHI B. 2020. A novel stochastic wastewater quality modeling based on fuzzy techniques. J Env Health Sci Eng 18(2): 1099-1120. doi:10.1007/s40201-020-00530-8.. Further, according to Hassani et al. (2019)HASSANI IE, MAZGUALDI CE & MASROUR T. 2019. Artificial Intelligence and Machine Learning to Predict and Improve Efficiency in Manufacturing Industry. arXiv URL http://arxiv.org/abs/1901.02256.
http://arxiv.org/abs/1901.02256 ...
hybrid Extremal-Micro genetic algorithm (Ex-µGA) used to tune the input parameters. Finally, three hybrid moving average-ANFIS (MA-ANFIS) models have been developed to validate the performance of proposed model. This will be used as a decision support tool to assist production managers in preventing various types of losses and then respond appropriately to the circumstance in order to maximize productivity. Figure 1 depicts the proposed system’s architecture, namely the Equipment Effectiveness Prediction System (EEPS).

Figure 1
Equipment effectiveness prediction system model.

ARCHITECTURE OF ANFIS MODEL DEVELOPMENT

ANFIS is exposed to be the best function approximator with fast convergence when compared to other hybrid neuro-fuzzy models (Akcayol 2004AKCAYOL MA. 2004. Application of adaptive neuro-fuzzy controller for SRM. Adv Eng Soft 35(3-4): 129-137. doi:10.1016/j.advengsoft.2004.03.005., Karaboga & Kaya 2019KARABOGA D & KAYA E. 2019. Training ANFIS by Using an Adaptive and Hybrid Artificial Bee Colony Algorithm (aABC) for the Identification of Nonlinear Static Systems. Arab J Sc Eng 44(4): 3531-3547. doi:10.1007/s13369-018-3562-y., Tiruneh et al. 2020TIRUNEH GG, FAYEK AR & SUMATI V. 2020. Neuro-fuzzy systems in construction engineering and management research. Autom Constr 119: 103348. doi:10.1016/j.autcon.2020.103348.). Mamdani’s fuzzy inference method is commonly used in developing prediction models in the Fuzzy inference system (FIS) analysis. However, when compared to Mamdani’s fuzzy inference method, Takagi-Sugeno (TS) model fuzzy inference achieves better results in model output interpretation (Wang & Chen 2008WANG WP & CHEN Z. 2008. A neuro-fuzzy based forecasting approach for rush order control applications. Expert Syst Appl 35(1-2): 223-234. doi:10.1016/j.eswa.2007.06.038.). Thus, TS fuzzy inference is used for model development in this study. The general architecture of the TS-based ANFIS model has been presented here for better comprehension. More details regarding ANFIS model can be found in many studies such as Chang (2008)CHANG BR. 2008. Resolving the forecasting problems of overshoot and volatility clustering using ANFIS coupling nonlinear heteroscedasticity with quantum tuning. Fuzzy Sets Syst 159(23): 3183-3200. doi:10.1016/j.fss.2008.04.003., Nguyen & Choi (2015)NGUYEN SD & CHOI SB. 2015. Design of a new adaptive neuro-fuzzy inference system based on a solution for clustering in a data potential field. Fuzzy Sets Syst 279: 64-86. doi:10.1016/j.fss.2015.02.012., Azimi et al. (2017)AZIMI H, BONAKDARI H, EBTEHAJ I, ASHRAF TALESH SH, MICHELSON DG & JAMALI A. 2017. Evolutionary Pareto optimization of an ANFIS network for modeling scour at pile groups in clear water condition. Fuzzy Sets Syst 319: 50-69. doi:10.1016/j.fss.2016.10.010.. ANFIS’ architecture is comprised of 5 fixed layers, each of which has a number of nodes specified as node function. This system is considered to have two inputs of ’x’ , ’y’ and an output of ’f’ to simplify the discussion. The two fuzzy if-then rules for a first order two rule TS type FIS can be computed as illustrated in equation 1-2.

If x is A 1 and y is B 1 then f 1 = a 1 x + b 1 y + c 1 = ω 1 (1)
If x is A 2 and y is B 2 then f 2 = a 2 x + b 2 y + c 2 = ω 2 (2)

In these, fuzzy sets are described as A1, A2, …, An. The output is described as ‘yi’ and output parameters mentioned as a,b,c. Furthermore, following are the FIS activities that are specified under the 5-layer architecture:

Layer 1: The first layer is referred to as the input layer. Each ‘i’ node in this layer is a square node. Equation 3 represents the node function. The membership function of ‘Ai’ is ‘O1i’ and it describes how much the provided ‘x’ satisfies the quantifier ‘Ai’.

O 1 i = μ Ai ( x ) = 1 1 + ( x c i a i ) 2 b i for i = 1 , 2 (3)

In this equation the output of ith node is represented as ‘O1i’. The inputs of ith node and membership function are represented as ‘x’, μAi(x) respectively. Likewise, the premise parameters are represented as ‘ai, bi, ci’.

Layer 2: This second layer consists of circular node and multiplication of input signals deliverers output, represented by ωi. The relation is expressed in equation 4.

O 2 i = ω i = μ Ai ( x ) μ Bi ( y ) for i = 1 , 2 (4)

Where, ωi – firing strength

Layer 3: This third layer is also comprised of circular nodes and is denoted by the symbol ω¯.

The ratio of the firing strength is calculated by the ith node using equations 5-6. This layer’s output is normalized firing strengths.

( ω ¯ ) = ( Σ ith node firing strength ) / ( Σ all rules firing strength ) (5)
O 3 i = ω i ¯ = ω 1 ω 1 + ω 2 for i = 1 , 2 (6)

Layer 4: In this layer, each interactive node is a square node that calculates the contributions of the ith node to overall output; the relationship is stated in equation 7.

O 4 i = ω i ¯ f i = ω i ¯ ( a i x + b i y + c i ) f o r i = 1 , 2 (7)

Layer 5: This layer’s only fixed node is a circular node indicated by notation ‘∑’. This calculates the network’s overall output by adding all of the rules’ contributions together. Equation 8 expresses the relationship.

O 5 i = Σ ω ¯ i f i = f = o v e r a l l o u t p u t for i = 1 , 2 (8)

The total output function defined as the linear combination of the consequent parameters in equation 9 with the values of the premise parameters constant.

f = ( ω 1 f 1 + ω 2 f 2 ) / ( ω 1 + ω 2 ) = ω ¯ f 1 + ω ¯ f 2 (9)

CONSTRUCTING AND VALIDATING OF HYBRID MA-ANFIS MODEL

The complete approach for the proposed model is presented in this section. The network models are developed and configured using the MATLAB 2018b ANFIS editor. After a trial and error method, the most appropriate network type is selected based on the features of the situation. The most appropriate training function for the problem is chosen as a result of a review of several training functions. ANFIS parameters have been optimized after numerous trials. Furthermore, the testing data has been utilized to evaluate the performance of prediction model, while training data is being used to train the network model. Finally, the proposed hybrid MA-ANFIS models are compared to determine which is the most effective.

Case selection

Human mobility is greatly influenced by the automobile industry. On-road and off-road vehicles are used to meet human transportation needs. In particular, buses are mostly used by the public sector for daily transportation. Bus body building is an essential milestone in the bus manufacturing process. The bus chassis is supplied by automobile manufacturers, while the body is built by vehicle bodybuilders. Medium-scale bus bodybuilding industries are back-born of the vehicle fabrication industry and this sector is strengthening by the massive number of labour forces. The Karur medium-scale bus bodybuilders play a vital role in bus bodybuilding sector of Tamil Nadu, India. A leading industry with three different cases was chosen to examine the feasibility of proposed research methodology. The selected medium-scale industry is specializing in fabricating, repairing, strip work, exterior and interior designing for bus manufacturing.

Case 1

The case-1 is used to the analysis of a real-world case involving of school buses body building. This first industrial case contributed a year’s worth of data, which was collected weekly from January to December 2019. There are total of 11 activities were observed to build the school bus. The fabrication of school bus included with more additional safety features. The initial measurement is done with actual equipment effectiveness. The initial measurement of mean OEE is 29 %, mean production loss due to stoppages (Lst) is 89 units, mean production loss due to reduced production speed (Lsp) is 564 units and mean quality loss (Lq) is 17 units. The losses are computed using equation 10-12.

L st = T LRS × S r (10)
L SP = ( T oaq ) ( S r × T o ) ( T oaq ) ( S p × T o ) (11)
L q = T ow T oaq (12)

Where TLRS represents realized stoppage duration, Sr represents realized production speed, Toaq represents total acceptable quantity of work units, To represents total operating time, Sp represents planned production speed and Tow represents total quantity of unit produced. In addition, total fabrication duration (DF) measured as 46 days, total fabrication cost (FC) is measured as Rs. 182029 and these are computed using equation 13-14.

D F = Q * P / D M H (13)

Where

F C = x = 1 n [ D LC x L u D F ] (14)

Where, ‘Q’ is quantity of work to be done, ‘DMH’ is daily man hours multiplied by daily straight hours (default: 8 hours), and ‘P’ is labour productivity. Likewise, ‘x’ represents activities carried out to complete fabrication, ‘n’ is number of activities, ‘DLC’ is direct labour cost, ‘Lu’ is labour usage.

Case 2

In the second case, the analysis of a real-world case involving of tour bus body building from the same industry selected in case-1. This industrial case also contributed a year’s worth of data, which was collected weekly with the same duration case-1. The same case-1 computation method was used to measure input values. The initial measurement is done with actual equipment effectiveness. The initial measurement of mean ‘OEE’ is 29 %, mean ‘Lst’ is 89 units, ‘Lsp’ is 564 units and ‘Lq’ is 17 units. Likewise, ‘DF’ is measured as 44 days, and ‘FC’ is measured as Rs. 179443.

Case 3

In the third case, the same industry as in case-1 was chosen. This case seems pertinent to the scenario of city bus body building. This industrial case also contributed weekly interval data and computation methods also follow same as mentioned in ‘case-1’. The initial measurement is done with actual equipment effectiveness. The initial measurement of mean ‘OEE’ is 44 %, mean ‘Lst’ is 89 units, ‘Lsp’ is 564 units and ‘Lq’ is 17 units. Likewise, ‘DF’ is measured as 41 days, and ‘FC’ is measured as Rs. 143780.

Input parameter selection

In most cases, standard production speed is taken into account while calculating OEE. This has the potential to conceal the true capability of the process. Similarly, this raises the OEE score artificially while concealing losses and limiting improvement. Furthermore, excessive data collecting by labours while performing task may result in work disruption. This could cause work to be slowed and productivity to be lost. On the other hand, using cutting-edge technologies such as IoT, for data collection while manufacturing, can help to solve these issues. This substitution technology is more expensive and suited to automated discrete manufacturing systems. However, this technology substitution is inapplicable when considering a manually operated or partially automated discrete manufacturing system of labour intensive medium scale bus body building works. Because technology substitution and labour management training require significant investment and raise direct labour costs. Alternatively, these complexities in OEE score improvement can be resolved by performing a critical analysis of production losses with respect to varying production speed in a real-world problem. This interesting combination occurs in the context of OEE, which has received little attention in the literature. Hence, the volatility of equipment effectiveness prediction using different loss functions is addressed in this study. Kuo & Lin (2010)KUO Y & LIN KP. 2010. Using neural network and decision tree for machine reliability prediction. Int J Adv Manuf Technol 50(9-12): 1243-1251. doi:10.1007/s00170-010-2593-7. predicted equipment effectiveness using only the loss function of equipment performance as a quantitative measure. Similarly, Engelmann et al. (2020)ENGELMANN B, SCHMITT S, MILLER E, BRÄUTIGAM V & SCHMITT J. 2020. Advances in machine learning detecting changeover processes in cyber physical production systems. J Manuf Mat Process 4(4): 108. doi:10.3390/jmmp4040108. used the loss function of equipment availability to predict equipment effectiveness. Hence, this critical analysis for selecting input parameter considers loss functions of availability, performance and quality with respect to varying production speed. Totally six input parameters that influence equipment effectiveness has been found as a result of these critical analyses. Machine availability (A), losses due to stoppages (Lst), machine performance (P), losses due to production speed (Lsp), machine quality (Q), and quality losses (Lq) were found in the literatures Kuo & Lin (2010)KUO Y & LIN KP. 2010. Using neural network and decision tree for machine reliability prediction. Int J Adv Manuf Technol 50(9-12): 1243-1251. doi:10.1007/s00170-010-2593-7., Ma et al. (2012)MA YH, ZHENG YB, ZHANG Y, JIANG Y & ZHAO LL. 2012. Diagnosis and prediction for equipment efficiency based on OEE in FCC M6 welding workshop. Appl Mech Mat 130-134: 1092-1096. doi:10.4028/www.scientific.net/AMM.130-134.1092. and Engelmann et al. (2020)ENGELMANN B, SCHMITT S, MILLER E, BRÄUTIGAM V & SCHMITT J. 2020. Advances in machine learning detecting changeover processes in cyber physical production systems. J Manuf Mat Process 4(4): 108. doi:10.3390/jmmp4040108.. Determining the probability (p-value) of obtaining a test statistic is essential for verifying the association between each input and output variable. The p-values are in the range of 0 to 1. The test statistic p-value must be at least as extreme as the computed value if the null hypothesis is true. In this study, the relationship between each input and output variable was verified using Pearson correlation analysis. This assessment has been based on the case-1 data sets. A commonly used alpha level of 0.01 was chosen for this study. The Pearson correlation matrix and p-values are illustrated in Table I. The results show that ‘Lst’, ‘A’, ‘P,’ and ‘Q’ are all positively related to ‘OEE.’ Similarly, the input parameters ‘Lsp’ and ‘Lq’ have an inverse relationship with the output parameter ‘OEE’. Correspondingly, the p-values of ‘Lst’, ‘Lsp’, and ‘Lq’ are 0, 0, 0 respectively. This confirms that the test statistic’s p-values are less than the chosen alpha level of 0.01. As a result, the null hypothesis is rejected. On the other hand, the p-values of ‘A’, ‘P’, and ‘Q’ have 0.382, 0.060, and 0.407, respectively. This statistic demonstrates that the p-values are greater than the alpha level of 0.01, implying that the null hypothesis is accepted. The results of the Pearson correlation analysis confirm the linear relationship between the input and output parameters. Furthermore, these results show that ‘Lst’, ‘Lsp’, and ‘Lq’ have low p-values. This indicates that the loss function variables ‘Lst’, ‘Lsp’, and ‘Lq’ are the best input parameters for OEE prediction.

Table I
Pearson correlation matrix with p-values..

Furthermore, the complexity of defining membership function also addressed in this research. The volatility of systematic method for defining the membership function is an important signal for predicting equipment effectiveness. Bekar et al. (2015)BEKAR ET, CAKMAKCI M & KAHRAMAN C. 2015. An ANFIS algorithm for forecasting overall equipment effectiveness parameter in total productive maintenance. J Mult-Valued Log Soft Comput 25(6): 535-554., used the triangle membership function to calculate membership degrees and developed an OEE prediction model. However, this solely considers triangle membership function. Hence, in this work, the triangle membership function, as well as two other membership functions from Gbell and Gauss, is used to compute membership degrees. In this study, three hybrid MA-ANFIS models have been developed to predict equipment effectiveness. Each model is made up of three independent input variables that lead to deterioration in equipment performance. The Case-1 model employs three loss functions and a triangle membership function to determine membership degrees. Case-2 uses the identical loss functions as case-1 and a Gbell-shaped membership function to calculate membership degrees. All given loss functions, as well as the Gauss membership function, are taken into account in the Case-3 model. Evaluation and error analysis are used to determine the optimal prediction model. As a result, various constraints for selecting maintenance operations have been proposed for assessing the decision-making of production managers.

Data sets selection

The designed system is initially assumed to perform OEE predictions on a weekly basis. The input variable is explained in full below. Recorded OEE in past ‘n’ periods, which including period t, t-1, t-2,... t-n, have been recognized as an important prediction signals for predicting equipment effectiveness in the next period t+1. The trend analysis is appropriate since the component deterioration trend is reported in a time serious form (Maurya et al. 2010MAURYA MR, PARITOSH PK, RENGASWAMY R & VENKATASUBRAMANIAN V. 2010. A framework for on-line trend extraction and fault diagnosis. Eng Appl Artif Intell 23(6): 950-960. doi:10.1016/j.engappai.2010.01.027.). In addition, Ma et al. (2012)MA YH, ZHENG YB, ZHANG Y, JIANG Y & ZHAO LL. 2012. Diagnosis and prediction for equipment efficiency based on OEE in FCC M6 welding workshop. Appl Mech Mat 130-134: 1092-1096. doi:10.4028/www.scientific.net/AMM.130-134.1092. also confirms that the data smoothing is a very useful technique to enhance data quality. Therefore, one of the realistic trend analysis models of moving average smoothing technique is used to smooth the input data. The moving average model uses an equation 15 for estimating the value.

y t = c + εt + θ 1 ε t 1 + θ 2 ε t 2 + + θ q ε t q (15)

Where, εt is noise, yt is weighted moving average forecast errors. The data has been further tuned to enhance the data quality. According to Hassani et al. (2019)HASSANI IE, MAZGUALDI CE & MASROUR T. 2019. Artificial Intelligence and Machine Learning to Predict and Improve Efficiency in Manufacturing Industry. arXiv URL http://arxiv.org/abs/1901.02256.
http://arxiv.org/abs/1901.02256 ...
, the optimization model of genetic algorithm used in this study. However, simple genetic algorithm probably produces an optimized data tuning, but, this will happen at a relatively slow rate. Hence, this is ineffective in real-world problem solving because it requires rapid convergence. Structure of the Micro Genetic Algorithm (μGA) is identical to that of simple GA, and it is used to work with small populations. Hence, it necessitates lesser function evaluations than traditional GA (Anderson et al. 2015ANDERSON BR, PRICE P, GUNAWIDJAJA R & EILERS H. 2015. Microgenetic optimization algorithm for optimal wavefront shaping. Appl Opt 54(6): 1485-1491. doi:10.1364/ao.54.001485.). The crossover operation is carried out in μGA. However, the mutation isn’t required, because sufficient variance had been maintained during micro population convergence. On the other hand, the self-organized critical idea is at the heart of the extremal optimization technique. The most unfavorable variables in a suboptimal response are gradually replaced by newly produced random variables in this method. The physical instinct to streamline also can serve as a source of inspiration for extremal optimization. This path is commonly followed by its forerunners, such as simulated annealing and GA (Majumdar et al. 2021MAJUMDAR A, DEBNATH T, BISWAS A, SOOD SK & BAISHNAB KL. 2021. An Energy Efficient e-Healthcare Framework Supported by Novel EO-μGA (Extremal Optimization Tuned Micro-Genetic Algorithm). Inf Syst Front 23(4): 1039-1056. doi:10.1007/s10796-020-10016-5.). In this study, the μGA is combined with the extremal optimization technique. The chromosome symbolizes the improving strategy, and the gene symbolizes the decision variable in this process. Therefore, the input parameter has been tuned using hybrid Extremal-micro genetic algorithm (Ex-µGA) to enhance the data quality. The Table II describes the data sets used for developing prediction models.

Table II
Descriptive data sets.

Training of fuzzy inference system parameters

The major purpose of creating this prediction model is to determine the effectiveness of equipment from the perspective of bus body builders. It’s worth noting that the impact of variables on forecasting isn’t always clear. The obtained data set, which consists of 52, 52, 52 data pairs, is being used to validate the neural system. The entire data set is split into two sections: training and testing. The training data sets represent approximately 80% of the data in the study, with the rest being utilized to test and improve the adaptive network, as evidenced in various studies (Hasanipanah et al. 2016HASANIPANAH M, JAHED ARMAGHANI D, KHAMESI H, BAKHSHANDEH AMNIEH H & GHORABA S. 2016. Several non-linear models in estimating air-overpressure resulting from mine blasting. Eng Comput 32(3): 441-455. doi:10.1007/s00366-015-0425-y., Mojtahedi et al. 2019MOJTAHEDI SFF, EBTEHAJ I, HASANIPANAH M, BONAKDARI H & AMNIEH HB. 2019. Proposing a novel hybrid intelligent model for the simulation of particle size distribution resulting from blasting. Eng Comput 35(1): 47-56. doi:10.1007/s00366-018-0582-x.). The testing data set, on the other hand, is used to determine whether the model has been over fit during training. One of the most beneficial TS inference systems is utilized to achieve better results and model output. When compared to Mamdani’s fuzzy inference, it has proven to be more accurate and easier to interpret (Wang & Chen 2008WANG WP & CHEN Z. 2008. A neuro-fuzzy based forecasting approach for rush order control applications. Expert Syst Appl 35(1-2): 223-234. doi:10.1016/j.eswa.2007.06.038., Ebtehaj & Bonakdari 2014EBTEHAJ I & BONAKDARI H. 2014. Performance Evaluation of Adaptive Neural Fuzzy Inference System for Sediment Transport in Sewers. Water Resour Manag 28(13): 4765-4779. doi:10.1007/s11269-014-0774-0.). The structure of Case-3 TS model is presented in Figure 2 as instance.

Figure 2
Structure of Takagi-Sugeno model.

Generally, subtractive clustering (SC), grid partitioning (GP) and fuzzy c-means (FCM) techniques are used to generate FIS structure. The ‘SC’ is fast as well as one-pass method of estimating the clusters in a set of data. This is useful when it’s unclear how many clusters should present for a given set of data (Chiu 1994CHIU SL. 1994. Fuzzy model identification based on cluster estimation. J Intell Fuzzy Syst 2: 267-278.). The GP divides the input space into various rectangle subspaces using a number of small fuzzy regions. The fuzzy rules, on the other hand, grow exponentially as the input variables increases (Moradi et al. 2019MORADI F, BONAKDARI H, KISI O, EBTEHAJ I, SHIRI J & GHARABAGHI B. 2019. Abutment scour depth modeling using neuro-fuzzy-embedded techniques. Mar Georesources Geotechnol 37(2): 190-200. doi:10.1080/1064119X.2017.1420113.). Furthermore, in FCM techniques, each data set is connected to a cluster to a degree specified by a membership grade. It demonstrates how to divide data points in a multidimensional space into 12 separate clusters (Bezdek 1981BEZDEK JC. 1981. Pattern Recognition with Fuzzy Objective Function Algorithms. Pattern Recognition with Fuzzy Objective Function Algorithms doi:10.1007/978-1-4757-0450-1., Safari et al. 2019SAFARI MJS, EBTEHAJ I, BONAKDARI H & ES-HAGHI MS. 2019. Sediment transport modeling in rigid boundary open channels using generalize structure of group method of data handling. J Hydrol 577: 123951. doi:10.1016/j.jhydrol.2019.123951., Lotfi et al. 2020LOTFI K, BONAKDARI H, EBTEHAJ I, DELATOLLA R, ZINATIZADEH AA & GHARABAGHI B. 2020. A novel stochastic wastewater quality modeling based on fuzzy techniques. J Env Health Sci Eng 18(2): 1099-1120. doi:10.1007/s40201-020-00530-8.). Therefore, this study used one of the most effective and appropriate grid partitioning techniques to generate membership functions.

Then, learning approach is one of the most important characteristics of neural networks. The hybrid learning function uses neural networks for training and learning to find fuzzy parameters. The fuzzy rules and membership functions (MF) can be fine-tuned through learning. The hybrid learning strategy for a neural fuzzy system consists of two phases: data-driven rule creation and back-propagation (BP) learning-based rule adjustment. Further, the back propagation and basic learning algorithm is combined with Least Square Estimates (LSE) to speed up the process. Back propagation is currently combined with the LSE in adaptive inference systems to update the parameters of membership functions. The hybrid learning methodology comprises a forward and backward run in each epoch. In the forward pass, the input data functions are used to compute each node response. The function will continue to run till the error value is computed. The error rates travelled from the output terminal to the input terminal during the backward pass, and the parameters were adjusted accordingly. The gradient descent methodology is now integrated with the BP-LSE technique to report the parameters of MFs. The input data functional procedure continues in the forward pass until the error measure evaluates each node output. In the backward pass, the error rates propagate and are updated using the gradient technique (Loganathan & Girija 2014LOGANATHAN C & GIRIJA KV. 2014. Investigations on Hybrid Learning in ANFIS. Int J Eng Res Appl 4(10): 31-37.). Hence, this study used LSE and BP gradient descent method.

In general, various membership function (MF) available in ANFIS editor, such as, trapezoidal, triangular, generalized bell, Gaussian curve, etc. Among this, three membership functions have been selected to develop proposed prediction model. The triangular membership function (Trimf) used to develop Case-1 model. This is the most primitive because of the triangular curve shapes. Three factors are used to identify three points: ‘a’ , ‘c’ for the feet, and ‘b’ for the tip of the curve (Raharja et al. 2021RAHARJA MA, DARMAWAN ID, NILAKUSUMAWATI DP & SUPRIANA IW. 2021. Analysis of membership function in implementation of adaptive neuro fuzzy inference system (ANFIS) method for inflation prediction. J Phys: Conf Ser 1722(1): 012005. doi:10.1088/1742-6596/1722/1/012005.). The generalized bell membership function (Gbellmf) used to develop case-2 model. The Gbellmf is a symmetrical shape that resembles a bell. This function has three parameters: ‘a’ specifies the breadth of the bell-shaped curve, ‘b’ is a positive integer, and ‘c’ determines the curve’s centre in the universe of discourse (Babanezhad et al. 2021BABANEZHAD M, BEHROYAN I, MARJANI A & SHIRAZIAN S. 2021. Artificial intelligence simulation of suspended sediment load with different membership functions of ANFIS. Neural Comput Appl 33(12): 6819-6833. doi:10.1007/s00521-020-05458-6.). The Gaussian membership function (Gaussmf) used to develop Case-3 model. This has been the most widely used MF for characterizing fuzzy systems due to its simple notation and smoothness. This does have its own set of benefits, including the fact that it is non-zero, smooth, and is defined by only two parameters that are optimized during the training process (Gholami et al. 2018GHOLAMI A, BONAKDARI H, EBTEHAJ I, GHARABAGHI B, KHODASHENAS SR, TALESH SHA & JAMALI A. 2018. A methodological approach of predicting threshold channel bank profile by multi-objective evolutionary optimization of ANFIS. Eng Geo 239: 298-309. doi:10.1016/j.enggeo.2018.03.030., Yaseen et al. 2017YASEEN ZM, EBTEHAJ I, BONAKDARI H, DEO RC, DANANDEH MEHR A, MOHTAR WHMW, DIOP L, EL-SHAFIE A & SINGH VP. 2017. Novel approach for streamflow forecasting using a hybrid ANFIS-FFA model. J Hydrol 554: 263-276. doi:10.1016/j.jhydrol.2017.09.007.). Table III depicts an overview of the training parameters for three selected cases.

Table III
Training parameters for ANFIS modeling.

RESULTS AND DISCUSSION

In the present work, three hybrid MA-ANFIS has been proposed to predict effectiveness of equipment. This is based on three independent losses and three type of (Tri, Gbell, and Gauss) membership function. Likewise, three statistical indices (RMSE, R2 and MAPE) are used to evaluate precisions of the predictive models.

Results of ANFIS parameters in relation to the proposed prediction models

Case 1

This proposed prediction model has been developed using triangular membership function. The initial FIS generated using Grid partitioning. The estimated value of each variable is used to define the discourse universe of each variable, which is then divided using Grid partitioning. The initial input membership function is [2 3 3] which generates 18 number of fuzzy rules. Then, triangular shape activation function computes fuzzy membership values and the function are characterized as equation 16.

y = f ( x ; a , b , c ) = max ( min ( x a b a , c x c b ) , 0 ) (16)

Where, ‘y’ is membership value, ‘x’ is input values for which to compute membership values. The ‘a, b, c’ are membership function parameter. In this ‘a, c’ define the feet and ‘b’ defines its peak of triangular membership function. These parameters can be set and updated by the hybrid learning algorithm during the training procedure. The quantity of membership functions for each input as well as the types of output membership functions followed constant functions. The proposed overall performance model’s ANFIS characteristics have 48 nodes under the best structure. Furthermore, there are a total of 39 pairs of parameters, with 18 linear parameters and 21 nonlinear parameters. In addition, 42 and 10 data samples have been used for training and testing, respectively. This proves that the network’s total number of data pairs is less than the number of training data pairs. The run time is 72 seconds. The ANFIS characteristic of the proposed model is presented in Table IV.

Case 2

This proposed prediction model has been developed using generalized bell membership function. The number and initial positions of cluster centers are estimated using Grid partitioning. The estimated value of each variable is used to define each variable’s discourse universe, which is then divided via Grid partitioning. The initial input membership function is [2 3 2] which generates 12 number of fuzzy rules. Then, bell shaped activation function is used to develop the ANFIS model in this study, and the function is characterized as equation 17.

y = f ( x ; r , s , t ) = 1 1 + | x t r | 2 s (17)

Where ‘y’ is membership value, ‘x’ is input values for which to compute membership values. The ‘r, s, t’ are the parameters that respectively control the slope, centre and width of the bell-shaped function. These parameters can be set and updated by the same as case 1. This learning methodology combines the LSM and the BP, followed by the gradient descent learning method. Initial bell shaped functions with specific parameters are assigned to each fuzzification neuron at the beginning. The neuronal function centre associated to input ‘xi’. Thus, domain of ‘xi’ is uniformly segregated. Further, the function widths and slopes are also chosen to allow enough overlap between the respective functions. The quantity of membership functions for each input as well as the types of output membership functions followed constant functions. The proposed overall performance model’s ANFIS characteristics have 44 nodes under the best structure, and 12 fuzzy rules have been constructed. Furthermore, there are a total of 33 pairs of parameters, with 12 linear parameters and 21 nonlinear parameters. This proves that the network’s total number of data pairs is less than the number of training data pairs. The run time is 76 seconds. The ANFIS characteristic of the proposed model is presented in Table IV.

Case 3

This proposed prediction model has been developed using Gaussian membership function. The number and initial positions of cluster centers are estimated using Grid partitioning. The minimum and maximum values of each variable are utilized to define each variable’s discourse universe, which is then divided via Grid partitioning. The initial input membership function is [2 2 3] which generates 12 number of fuzzy rules. Then, Gaussian function is used to develop the ANFIS model in this study, and the function is characterized as equation 18.

y = f ( x ; σ , c ) = e ( x c ) 2 2 σ 2 (18)

Where ‘y’ is membership value, ‘x’ is input values for which to compute membership values. The ‘σ’ is standard deviation, ‘c’ is mean. These parameters can be set and updated by the same as case 1 learning methodology. The quantity of membership functions for each input as well as the types of output membership functions followed constant functions. The proposed overall performance model’s ANFIS characteristics have 44 nodes under the best structure, and 12 fuzzy rules have been constructed. Furthermore, there are a total of 26 pairs of parameters, with 12 linear parameters and 14 nonlinear parameters. This proves that the network’s total number of data pairs is less than the number of training data pairs. The run time is 70 seconds. The ANFIS characteristic of the proposed model is presented in Table IV.

Table IV
Characteristics of best structure model and corresponding ANFIS information.

Discussion of proposed MA-ANFIS prediction models evaluation

The trained fuzzy inference system then concern to testing data with various parameter combinations to evaluate the neural system’s generalization capabilities and avoid over fitting the training data set. Based on the testing error of each combination of parameters, the optimum configuration for the chosen cases is determined. The variables of proposed model are assumed as constant function for several past ‘p’ measured observations, ‘q’ prediction observation and ‘n’ number of samples. The performances among these ANFIS models are statistically assessed by root mean squared error (RMSE), the coefficient of determination (R2) and mean absolute percentage error (MAPE). The relation is represented in equation 19-21.

Root mean squared error = i = 1 n ( p i q i ) 2 n (19)
R squared = 1 i = 1 n ( p i q i ) ¯ 2 i = 1 n ( p i p i ) ¯ 2 (20)
M A P E = 1 n i = 1 n | p i q i p i * 100 | (21)

The OEE measured and OEE predicted values are statically evaluated by means of error analysis. This is worth noting that a prediction with an RMSE of 0; R2 of 1; and MAPE of 0 is a perfect model. Table V shows the error analysis for the selected three cases predictive models using training and testing datasets. Based on error analysis the case 3 model possess minimal RMSE of 0.87, high predictive accuracy of 98.9 percent and minimal MAPE of 1.5 percent. Thus, the case 3 predictive model (Gaussianmf based model) has the best performance for the OEE estimation. Further, the Figure 3 a-c, shows the values of OEE measured versus OEE predicted by the predictive models.

Figure 3
Testing data- Measured versus predicted R2 accuracy comparison.
Table V
Error analysis comparison.

As shown in these figures the prediction through hybrid MA-ANFIS model case 3 are very accurate and closer to measured values. The high R2 values shown by case 3 developed models as compared to other models indicate better prediction capability. The optimal values of the preceding and subsequent parameters for the Case-3 model are shown in Table VI.

Table VI
Optimal parameters of Gaussian membership functions.

Comparison with other existing work

The performances of proposed models are compared with existing prediction models to validate the prediction ability. Three previous models were chosen to compare the current proposed models. Existing models simply employed one statistical index to validate results. However, the findings of this investigation have been validated using three statistical indices. The Table VII shows the comparison of proposed and existing models test results with respect to RMSE, R2 and MAPE. These findings of error measures confirms that the accuracy of Case-3 proposed model is outperform than other models in predicting OEE. In combination with existing prediction research, this study demonstrates that the Gaussian membership function based hybrid MA-ANFIS model is capable of handling production losses for accurate equipment effectiveness prediction.

Table VII
Comparison with other existing work.

Statistical hypothesis test

This statistical test provides a mechanism for making quantitative decisions about the prediction process. The superiority case-3 prediction model results are selected for doing statistical hypothesis test. The intent is to determine whether there is enough evidence to accept a hypothesis about the prediction process. The paired t-test is selected for conducting hypothesis test, because the mean difference between OEE measured and OEE predicted are paired observations. The paired t-test computes the difference between each before and after pair of measurements. This calculates the mean of these OEE changes, and reports whether this mean of the difference is statistically significant. The test uses the 95 % confidence interval (CI) and following hypothesis:

H0: µ (OEE measured) = µ (OEEpredicted) (OEE prediction value is closer to measured value)

H1: µ (OEE measured) ≠ µ (OEEpredicted) (OEE prediction value is not closer to measured value)

The Table VIII demonstrates the paired T-test results. The calculated probability of obtaining the observed sample data is 0.960. This confirms that accept the null hypothesis, which means the measured and predicted OEE values are closer. This determined quantitative statistical evidence provides decisions to accept the hypothesis of the prediction process.

Table VIII
Paired T-test results.

Practical implication of this research

The greatest analysis of prediction model is to extracting valuable knowledge from raw data and assists production managers to determine the cut-off percentage of equipment effectiveness for selecting appropriate maintenance technique. The average equipment effectiveness value predicted by the superior Gaussianmf based hybrid MA-ANFIS model (case-3) is 34 percentages. However, this OEE predicted value is lower than the fair (65%) and world class (85%) OEE values (OEE predicted > OEEfair > OEEworld class) and this is required of further improvement. Hence, this work initiates failure analysis as well as countermeasures to reduce production losses and to enhance OEE values. The subsystem’s failure patterns aid production managers in understanding failure behavior and designing preventative maintenance programs. Finally this work proposes cut-off values of OEE for selecting appropriate maintenance technique. This is also provides higher accuracy in the evaluation of losses and prediction of effectiveness in a complex environment. Initially, losses due to failures and causes of machine centers have been investigated from a human, electrical, and mechanical perspective. The Table IX showcases the investigation statistics.

Table IX
Failure and repair statistics.

According to data, human failures contribute to the high number of failures in the machine centre. This kind of failures is caused by commission flaws such as a lack of attention, inadequate tool and material arrangements, and so on. Human failures have a mean time to repair (MTTR) of 0.18 hour and a mean time between failures (MTBF) of 1.6 hours. Next to human failure, the mechanical failures have been analyzed. Mechanical failures accounted for 18 percent of all failures, with one occurring every three hours on average and a 0.35-hour average repair time. Ultimately, mechanical failures are linked to failures of electrode holders, circuits, wires, liners, hoses, regulators, and inside guns, resulting in a stoppage. Electrical failures have also been studied, in addition to mechanical failures. Electrical problems accounted for 48 percent of all failures. Electrical failure had the second-shortest MTBF, with one failure per 1.13 hours and the shortest repair time of 0.11 hours. These evaluations are essential to analyze the current situation in order to enforce the appropriate counter measures. The first line of counter measure is self-maintenance procedures like cleaning, rod storage, and equipment management. This means maintaining the right equipment by following the principles of the Japanese 5’S. To regularly update the machine, preventive maintenance (PM) has been encouraged. The electrical maintenance countermeasures include weekly power source monitoring by blowing air into the unit, which reduces overheating and intermittent arc output. Problems such as excessive circuit resistance would also be caused by the cable and electrode holder, so the connections, circuits, wires, liners, hoses, regulators and inside guns should be regularly inspected. Direct replacement ones with instant excess wear. Of course, with blended coolant from the manufacturer, the coolant level is regulated to solve the problem related to sludge build-up. Concentrate primarily on generator welding engines such as oil change, filter, air cleaner and filter for fuel. Lastly, conduct the load test to ensure that the welding output is precise. A failure monitoring template has been developed with a view to supporting maintenance activities. The model consists of a critical complete record of often missing details related to relevant specifics, such as machine number, the period of failure detection, mean repair time and the mean time between repair failures. This lower buffer inventories result in low production losses and shorter lead times because tasks do not have to wait as long in lines. So, the company’s performance and adaptability have improved, making it more competitive. Increases in quality entail reduced scrap and rework, resulting in a higher quality rate. In addition, the countermeasures have been estimated for each selected parameters and OEE values are estimated using proposed prediction model. Likewise, the cut-off value of OEE is proposed as follows,

OEE cut-off value Equipment status Proposed maintenance strategy
85 ≤ OEE predicted Excellent Follow same PM technique
65 ≤ OEE predicted ≤ 85 Good Follow same PM technique
55 ≤ OEE predicted ≤ 65 Normal Implement proposed PM technique
OEE predicted ≤ 55 Fail Rigorous implementation of proposed PM method

The OEE prediction using proposed model are described in Figure 4. The predicted average OEE value is 66 percentages. As a result, the OEE has risen from 34 to 66 percent. Thus, the equipment effectiveness reached fair OEE of 65%. Now the OEE value is higher than fair value and lower than world class value (OEE fair > OEEpredicted > OEEworld class). In order to reach world-class performance of 85 % OEE, the production managers must continue these countermeasures effectively. This prediction can assist production managers to choose appropriate maintenance strategy especially for manual and semi automated production environment scenario.

Figure 4
OEE prediction using proposed model.

Cost effectiveness of proposed study

The superior proposed model (case-3) is used to analyze a real world case which is fabrication of city buses. The production manager could be able to enhance labour productivity by using the integrated improvement method. The results show that when improvement factor (EEI) performs the maximum output valued at average of 86 percentages. Due to the improvement factor the average optimal LTTC * is reduced from 0.56 to 0.46 (man hours required to complete unit quantity of work). The optimal labour productivity (LP *) increased to 2.84 from 1.95 units per man hour. The total fabrication duration (DF) could be reduced from 41 days to 27 days. The improvement cost of this work is calculated as 19,000 rupees. The optimal direct labour cost associated with reduced duration is 120486. The methods surveyed are compared in terms of OEE, productivity, time, and cost effectiveness in Figure 5a-d.

Figure 5
Effectiveness of the proposed methods surveyed.

CONCLUSION

There are enormous challenges in today’s volatile industrial environment to improve operational efficiency in order to boost competitiveness and overall business success. In the dynamic environment, medium-scale bus body builders are struggling to retain existing customers and gain new customers. Bus body builders who produce high-quality vehicles will be able to thrive in core manufacturing. These builders work hard to increase their competitive advantages by focusing on the efficiency of their equipment. On the other hand, Precision OEE prediction is critical in these kinds of extremely difficult manufacturing environments. This study combines a moving average (MA) model and hybrid Extremal- micro genetic algorithm (Ex-µGA) model with an adaptive neural fuzzy inference system (ANFIS) to develop three different prediction models for equipment effectiveness prediction. Various stages of evaluation were carried out in order to determine the best ANFIS parameters and validate the ability of the hybrid MA-ANFIS model. The empirical results clearly show that the performance of the Gaussianmf-based hybrid MA-ANFIS model (case-3) outperforms other developed models. The proposed system assists production managers by providing decision support for quantifying and predicting OEE performance. The proposed equipment effectiveness prediction can aid in the effective implementation of maintenance management programs in the manufacturing environment. The results of implementing the superiority model in a bus body building firm show significant improvements in labour productivity, machine performance, fabrication time, and fabrication cost. According to the real-world analysis, bus body building production managers must track failures on a regular basis. Otherwise, failures will have increased significantly in comparison to the projected output.

Limitation and future direction of proposed study

This section discusses the limitations of the current study as well as an attractive future direction. The present hybrid MA-ANFIS superiority model based on Gaussianmf has the ability to immediately improve task performance. However, it has only used grid partitioning clustering and the Ex-µGA as input parameter tuning technique. Hence, few other clustering techniques and input parameter tuning methods can be given more strength in OEE prediction. In addition, the production managers have to predict some other flexible preventive maintenance technique to accelerate equipment performance within the budgeted range.

  • AHUJA IPS & KHAMBA JS. 2008. Total productive maintenance: literature review and directions. Int J Qual Reliab Manag 25(7): 709-756.
  • AKCAYOL MA. 2004. Application of adaptive neuro-fuzzy controller for SRM. Adv Eng Soft 35(3-4): 129-137. doi:10.1016/j.advengsoft.2004.03.005.
  • ANDERSON BR, PRICE P, GUNAWIDJAJA R & EILERS H. 2015. Microgenetic optimization algorithm for optimal wavefront shaping. Appl Opt 54(6): 1485-1491. doi:10.1364/ao.54.001485.
  • ANVARI F, EDWARDS R & STARR A. 2011. Evaluation of overall equipment effectiveness based on market. J Qual Maint Eng 16(3): 256-270.
  • AZIMI H, BONAKDARI H, EBTEHAJ I, ASHRAF TALESH SH, MICHELSON DG & JAMALI A. 2017. Evolutionary Pareto optimization of an ANFIS network for modeling scour at pile groups in clear water condition. Fuzzy Sets Syst 319: 50-69. doi:10.1016/j.fss.2016.10.010.
  • BABANEZHAD M, BEHROYAN I, MARJANI A & SHIRAZIAN S. 2021. Artificial intelligence simulation of suspended sediment load with different membership functions of ANFIS. Neural Comput Appl 33(12): 6819-6833. doi:10.1007/s00521-020-05458-6.
  • BEKAR ET, CAKMAKCI M & KAHRAMAN C. 2015. An ANFIS algorithm for forecasting overall equipment effectiveness parameter in total productive maintenance. J Mult-Valued Log Soft Comput 25(6): 535-554.
  • BEZDEK JC. 1981. Pattern Recognition with Fuzzy Objective Function Algorithms. Pattern Recognition with Fuzzy Objective Function Algorithms doi:10.1007/978-1-4757-0450-1.
  • CHANG BR. 2008. Resolving the forecasting problems of overshoot and volatility clustering using ANFIS coupling nonlinear heteroscedasticity with quantum tuning. Fuzzy Sets Syst 159(23): 3183-3200. doi:10.1016/j.fss.2008.04.003.
  • CHIU SL. 1994. Fuzzy model identification based on cluster estimation. J Intell Fuzzy Syst 2: 267-278.
  • CHOU TC, HSU LL, YEH YJ & HO CT. 2005. Towards a framework of the performance evaluation of SMEs’ industry portals. Ind Manag Data Syst 105(4): 527-544. doi:10.1108/02635570510592406.
  • EBTEHAJ I & BONAKDARI H. 2014. Performance Evaluation of Adaptive Neural Fuzzy Inference System for Sediment Transport in Sewers. Water Resour Manag 28(13): 4765-4779. doi:10.1007/s11269-014-0774-0.
  • ENGELMANN B, SCHMITT S, MILLER E, BRÄUTIGAM V & SCHMITT J. 2020. Advances in machine learning detecting changeover processes in cyber physical production systems. J Manuf Mat Process 4(4): 108. doi:10.3390/jmmp4040108.
  • FATTAH J, EZZINE L & LACHHAB A. 2017. Evaluating the performance of a production line by the overall equipment effectiveness: An approach based on best maintenance practices. Int J Eng Res Afr 30: 181-189. doi:10.4028/www.scientific.net/JERA.30.181.
  • GHOLAMI A, BONAKDARI H, EBTEHAJ I, GHARABAGHI B, KHODASHENAS SR, TALESH SHA & JAMALI A. 2018. A methodological approach of predicting threshold channel bank profile by multi-objective evolutionary optimization of ANFIS. Eng Geo 239: 298-309. doi:10.1016/j.enggeo.2018.03.030.
  • HASANIPANAH M, JAHED ARMAGHANI D, KHAMESI H, BAKHSHANDEH AMNIEH H & GHORABA S. 2016. Several non-linear models in estimating air-overpressure resulting from mine blasting. Eng Comput 32(3): 441-455. doi:10.1007/s00366-015-0425-y.
  • HASSANI IE, MAZGUALDI CE & MASROUR T. 2019. Artificial Intelligence and Machine Learning to Predict and Improve Efficiency in Manufacturing Industry. arXiv URL http://arxiv.org/abs/1901.02256.
    » http://arxiv.org/abs/1901.02256
  • KARABOGA D & KAYA E. 2019. Training ANFIS by Using an Adaptive and Hybrid Artificial Bee Colony Algorithm (aABC) for the Identification of Nonlinear Static Systems. Arab J Sc Eng 44(4): 3531-3547. doi:10.1007/s13369-018-3562-y.
  • KUO Y & LIN KP. 2010. Using neural network and decision tree for machine reliability prediction. Int J Adv Manuf Technol 50(9-12): 1243-1251. doi:10.1007/s00170-010-2593-7.
  • KURSCHEIDT NETTO RJ, SANTOS EAP, DE FREITAS ROCHA LOURES E & PIEREZAN R. 2017. Using Overall Equipment Effectiveness (OEE) to Predict Shutdown Maintenance. In: AMORIM M, FERREIRA C, VIEIRA JR M & PRADO C (Eds), Engineering Systems and Networks. Lecture Notes in Management and Industrial Engineering. Berlin: Springer, p. 13-21. doi:10.1007/978-3-319-45748-2_2.
  • LJUNGBERG Õ. 1988. Measurement of overall equipment effectiveness as a basis for TPM activities. Int J Oper & Prod Manag 18(37): 510-528. URL http://dx.doi.org/10.1108/01443579810206334%0.
  • LOGANATHAN C & GIRIJA KV. 2014. Investigations on Hybrid Learning in ANFIS. Int J Eng Res Appl 4(10): 31-37.
  • LOTFI K, BONAKDARI H, EBTEHAJ I, DELATOLLA R, ZINATIZADEH AA & GHARABAGHI B. 2020. A novel stochastic wastewater quality modeling based on fuzzy techniques. J Env Health Sci Eng 18(2): 1099-1120. doi:10.1007/s40201-020-00530-8.
  • MA YH, ZHENG YB, ZHANG Y, JIANG Y & ZHAO LL. 2012. Diagnosis and prediction for equipment efficiency based on OEE in FCC M6 welding workshop. Appl Mech Mat 130-134: 1092-1096. doi:10.4028/www.scientific.net/AMM.130-134.1092.
  • MAJUMDAR A, DEBNATH T, BISWAS A, SOOD SK & BAISHNAB KL. 2021. An Energy Efficient e-Healthcare Framework Supported by Novel EO-μGA (Extremal Optimization Tuned Micro-Genetic Algorithm). Inf Syst Front 23(4): 1039-1056. doi:10.1007/s10796-020-10016-5.
  • MAURYA MR, PARITOSH PK, RENGASWAMY R & VENKATASUBRAMANIAN V. 2010. A framework for on-line trend extraction and fault diagnosis. Eng Appl Artif Intell 23(6): 950-960. doi:10.1016/j.engappai.2010.01.027.
  • MOJTAHEDI SFF, EBTEHAJ I, HASANIPANAH M, BONAKDARI H & AMNIEH HB. 2019. Proposing a novel hybrid intelligent model for the simulation of particle size distribution resulting from blasting. Eng Comput 35(1): 47-56. doi:10.1007/s00366-018-0582-x.
  • MORADI F, BONAKDARI H, KISI O, EBTEHAJ I, SHIRI J & GHARABAGHI B. 2019. Abutment scour depth modeling using neuro-fuzzy-embedded techniques. Mar Georesources Geotechnol 37(2): 190-200. doi:10.1080/1064119X.2017.1420113.
  • NAKAJIMA S. 1988. Introduction to total productive maintenance. Student Study Guide. New York: Productivity Press, p. 129.
  • NGUYEN SD & CHOI SB. 2015. Design of a new adaptive neuro-fuzzy inference system based on a solution for clustering in a data potential field. Fuzzy Sets Syst 279: 64-86. doi:10.1016/j.fss.2015.02.012.
  • NIKOLOPOULOS K, METAXIOTIS K, LEKATIS N & ASSIMAKOPOULOS V. 2003. Integrating industrial maintenance strategy into ERP. Indust Manag Data Syst 103(3-4): 184-191. doi:10.1108/02635570310465661.
  • OK SC & SINHA SK. 2006. Construction equipment productivity estimation using artificial neural network model. Construction Management and Economics 24(10): 1029-1044. doi:10.1080/01446190600851033.
  • PRASETYO YT & VEROYA FC. 2020. An Application of Overall Equipment Effectiveness (OEE) for Minimizing the Bottleneck Process in Semiconductor Industry. In: 2020 IEEE 7th International Conference on Industrial Engineering and Applications, ICIEA 2020, p. 345-349. doi:10.1109/ICIEA49774.2020.9101925.
  • RAHARJA MA, DARMAWAN ID, NILAKUSUMAWATI DP & SUPRIANA IW. 2021. Analysis of membership function in implementation of adaptive neuro fuzzy inference system (ANFIS) method for inflation prediction. J Phys: Conf Ser 1722(1): 012005. doi:10.1088/1742-6596/1722/1/012005.
  • RELKAR AS & NANDURKAR KN. 2012. Optimizing & analysing overall equipment effectiveness (OEE) through design of experiments (DOE). Procedia Eng 38: 2973-2980. doi:10.1016/j.proeng.2012.06.347.
  • SAFARI MJS, EBTEHAJ I, BONAKDARI H & ES-HAGHI MS. 2019. Sediment transport modeling in rigid boundary open channels using generalize structure of group method of data handling. J Hydrol 577: 123951. doi:10.1016/j.jhydrol.2019.123951.
  • SHARMA R. 2019. Overall equipment effectiveness measurement of TPM manager model machines in flexible manufacturing environment: a case study of automobile sector. Int J Prod Qual Manag 26(2): 206-222. doi:10.1504/ijpqm.2019.10018982.
  • TIRUNEH GG, FAYEK AR & SUMATI V. 2020. Neuro-fuzzy systems in construction engineering and management research. Autom Constr 119: 103348. doi:10.1016/j.autcon.2020.103348.
  • WANG FK & LEE W. 2001. Learning curve analysis in total productive maintenance. Omega 29(6): 491-499. doi:10.1016/S0305-0483(01)00039-1.
  • WANG WP & CHEN Z. 2008. A neuro-fuzzy based forecasting approach for rush order control applications. Expert Syst Appl 35(1-2): 223-234. doi:10.1016/j.eswa.2007.06.038.
  • YASEEN ZM, EBTEHAJ I, BONAKDARI H, DEO RC, DANANDEH MEHR A, MOHTAR WHMW, DIOP L, EL-SHAFIE A & SINGH VP. 2017. Novel approach for streamflow forecasting using a hybrid ANFIS-FFA model. J Hydrol 554: 263-276. doi:10.1016/j.jhydrol.2017.09.007.
  • YUNIAWAN D, ITO T & MOHAMAD ED. 2013. Measurement of Production Line Performance Based on OEE Extension. In: The Proceedings of Manufacturing Systems Division Conference, 57-58 p. doi:10.1299/jsmemsd.2014.57.

Publication Dates

  • Publication in this collection
    12 Dec 2022
  • Date of issue
    2022

History

  • Received
    12 Apr 2021
  • Accepted
    28 Oct 2021
Academia Brasileira de Ciências Rua Anfilófio de Carvalho, 29, 3º andar, 20030-060 Rio de Janeiro RJ Brasil, Tel: +55 21 3907-8100 - Rio de Janeiro - RJ - Brazil
E-mail: aabc@abc.org.br